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In this work, we perform a comparison of the employ of distinct gluon propagators 
with the experimental data in diffractive processes, pp elastic scattering and light me- 
son photo-production. The gluon propagators are calculated through non-perturbative 
methods, being justified their use in this class of events, due to the smallness of the mo- 
mentum transfer. Our results are not able to select the best choice for the modified gluon 
propagator among the analyzed ones, showing that the application of this procedure in 
this class of high energy processes, although giving a reasonable fit to the experimental 
data, should be taken with same caution. 

show of the results for gluon propagator can perform with care, . 
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1. Introduction 

The challenging description of the diffractive processes in high energy collisions 
in the framework of Quantum Chromodynamics has been requiring the attention 
over the last three decadei P^'H The main object of study is the Pomeron, which 
was originally proposed in Regge theory to explain the hadron collisions in the 
high energy limit. With the advent of QCD, the Pomeron was described by the 
exchange of a pair of gluons, in a first approximation or a ladder of gluons, the 
BFKL Pomerorffl 

However, in the kinematic region of small momentum transfer (called infrared 
region), we have divergences, due to the Landau poles in the gluon propagator. 
To remedy this problem, Landshoff and NachtmanrPI proposed a model for the 
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Pomeron, in which the perturbative gluon propagator is replaced by a modified one 
in the infrared region by the inclusion of non-perturbative effects. Following the 
years, this model (or at least the main idea) was employed in several applications 
in diffractive scattering, as, for example, pp elastic scattering, light meson photo- 
production. 

As we review in the next section, there are several solutions for the gluon prop- 
agator in the literature, each one with a particular choice of parameters (specially 
the massive ones). In this work, we choose to employ mainly the solutions which 
have a dynamic gluon mass, which can have, or not, a finite mass for the gluon 
when the momentum is zero, but matches the perturbative behavior in large mo- 
mentum transfer. There are also a direct connection of the gluon propagator"^ and 
the frozen of the running coupling constant, which means a finite value of a s for 
null momentum scale. 

The leitmotiv of this work is to compare the different forms of the gluon propa- 
gator in an independent framework: the diffractive processes. There is a possibility 
to perform just a fit to the data, adjusting the massive parameters, however the 
idea here is to compare the original values found for the parameters, to make a true 
comparison among the propagators. This approach is supported by the previous 
attempts in this direction (see, for example, Refs.EHEllZIEIand references therein). 

In the section El we briefly review the methods to obtain the modified gluon 
propagators and presents the solutions with a dynamic gluon mass. The models 
to describe the diffractive scattering collisions are presented in the section [3] and 
their results are shown and discussed in the section We present our conclusions 
in section |SJ 

2. The Gluon propagator and the running coupling constant 

There are several methods to obtain the gluon propagator by using non-perturbative 
methods. The more popular are the Dyson-Schwinger equations and the lattice field 
theory. 

The DSE's are an infinite system of non-linear coupled integral equations which 
relate the different Green functions of a quantum field theory. In the case of QCD, 
due to the complexity of the system to be solved, there are several analytical ap- 
proximations in the literature to obtain a solution. A discussion among the different 
methods and approximations employed can be found in Refs. 03 1101 

The lattice field theory, in opposition to the DSE, is a powerful numerical 
method, where the main idea is that the divergences of the theory are regular- 
ized by the discretization of the space-time. Otherwise, the lattice approximation 
is not free of problems: finite size of the lattice, spacing between the sites of the 
lattice and fermion simulation. 

Due to the approximations employed in the resolution of the problems pointed 
above, there are many different solutions for both methods in the literature. We 
focus here on the solutions which have the particular property of a dynamical gluon 
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mass: a gluon p ropaga tor with this feature was used in successful descriptions of 
several processegSEEl The first solution for the gluon propagator which will be 
used in this work is the Cornwall's solution^! from DSE's in the axial gauge with 
the use of a technique of re-summation of Feynman diagrams, 



with 



D c \q 2 )= [q 2 +m^(q 2 )] bg 2 In 



m^q 2 ) = 



q 2 + 4m 2 7 (q 2 ) 



A 2 




-i -12/11 



(1) 



(2) 



where b = 33/(487r 2 ) is the leading order coefficient of the /3 function and g = 
47ra s (0), with a s the running coupling constant. The value for m g lies i n th e interval 
500 ± 200 MeV for A QCD = 300 M eV giv en by the method of solution^! Another 
solution is obtained by Habel et a/™D31 using a different type of approximation, 
which gives the following propagator, 

where b is a parameter (for its value, see below). Gorbar and Natale^l use the 
propagator 



[^^(q 2 )] 1 = q 2 + A*2© W - q 2 ) + 4 Q (q 2 - W 



(4) 



to calculate the vacuum QCD energy through effective potentials, where O is the 
step function, fi g = 0.61149 GeV is related with the gluon condensate and £' = 
0.9666797 is a calculated parameter. A recent result for the gluon propagator, ob- 
tained from the DSE's in the Mandelstam approximation, is due to Aguilar and 
Natale^, is given by the expression 

D AN {q 2 ) = 2 ^ * 2V M 2 (q 2 ) = ^^, (5) 

where to 2 , = 0.99 GeV 2 for the QCD scale A QCD = 335 MeV. 

Another recent result is the so lution of DSE's in the Landau gauge calculated 
by Alkofer and collaborator J1U | 17 | 18 | 19 | -^^3 gl uon propagator is given 



°AL{q 2 ) = -^ 



(«i AL) (q 2 ) 



(6) 



^qcD + q- 

where u = 2.5, A QCD = 510 MeV, k « 0.595, 7 = -13/22 and ai AL) (q 2 ) is the 
running coupling constant, 



ai AL )(q 2 ) 
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l + (q 2 /A 2 QCD 
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,(7) 



2, 2008 6:34 WSPC/INSTRUCTION FILE ws-newpgm 



4 M. B. Gay Ducati, Werner K. Sauter 




w* io" 4 io" 2 10° io 2 10"' io° io' 10 

k 2 (GeV 2 ) k 2 (GeV 2 ) 



(a) (b) 

Fig. 1. (a) Distinct modified gluon propagators compared with the pcrturbative one. The prop- 
agators are indicated in the plot legend, (b) Comparison between different results for the running 
coupling constant. In the detail, two results of the Cornwall frozen coupling with two different £'s. 
The parameters of the curves are given in the text. 



where a s (0) « 2.972 and fa = 11. 

The momentum dependences of the gluon propagators are displayed in the figure 
Q in comparison with the perturbative one. The most remarkable feature is the 
infrared behavior (k — > 0), where exists a finite value in the case of the Cornwall 
propagator, a zero value, as in the Habel/Gribov propagator, and the divergent 
perturbative propagator. Another remarkable fact is the same ultraviolet behavior 
(k — > oo) of the modified propagators and the perturbative one. 

Nevertheless, there is a narrow relation between the gluon propagator and 
the running coupling, widely discussed in Refs. |SJ The phenomenological ap- 
plications o f a fr ozen coupling constant include heavy quarkonia decays^", meson 
form factorJ^El an d other observableJ^. In some of the above results (specially in 
Ref. Hip, the data available are described when the modified gluon propagator and 
the frozen coupling constant are used at the same time. 

The analytical expre ssion for the running coupling is obtained following the 
Cornwall's solution I I namely, 

c(^) = 4,/ ftl .( ^±Mg) ), (8) 

where the explicit momentum dependence is the same as the one loop perturbative 
calculation, but with an extra term, $mp(fe 2 ), resulting in the frozen infrared be- 
havior where £ is a parameter that varies with the p ropagator and measures the 
frozenness of the coupling. In the original woril ' \ the value of this parameter 
is obtained as £ w 4 or larger. The massive term, m-p(k 2 ), also depends on the 
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gluon propagator employed, being the dynamical mass term of the propagator. For 
example, in the case of the Aguilar and Natale propagator, eq. this term is 
m-p{k 2 ) = A4 2 (k 2 ). The exception is the Alkofer case, with a proper expression for 
the running coupling constant, eq. (JJJ. 

In the fig.n we present the different a s 's (given by eq.JSJ), where the freezing 
of the running coupling constant, the matching of the frozen running coupling, and 
the behavior of the perturbative result in high momentum are very clear. The last 
result comes from the behavior of the massive term (m-p(fc 2 )) which goes to zero 
when the momentum is large in all the propagators employed in this work. The 
values of the propagator parameters are from the original works or assumed in the 
region of validity of the models, as in case of the Cornwall propagator. The figure 
also shows the £ parameter dependence of the Cornwall solution, eq. Q . As seeing 
in fig. n there is a small difference for distinct £ values. 



3. Diffractive scattering processes 

The framework in which the different results for the gluon propagator are compared 
is the diffractive scattering. This class of QCD processes is characterized by the 
exchange of the Pomeron and the consequent rapidity gaps in the datcPl In the 
regime of small momentum transfer, it is a good laboratory to test the modified 
gluon propagators, due to the range of validity of the approximation employed in 
the derivation of these objects. 

The most employed model in this case is the Landshoof-Nachtmann on where 
the Pomeron is modeled by the exchange of two non-perturbative gluons (with 
modified gluon propagators). In the original version of the LN model, the parameters 
are fixed relating the Pomeron strength radius and the Pomeron effective mass. A 
disadvantage of the LN model is the energy independent cross sections, due to the 
two gluon approximation. 

We will analyze two processes in which the small momentum contributions are 
relevant and are mediated by the Pomeron exchange: the pp scattering and the 
electro-production of light mesons; althought already present in the literature, not 
with the emphasis assumed in this present work. 

The model employed to describe t he p roton-proton elastic scattering in this 
work was proposed by Cudell and Rosa^, based on the LN two-gluon model for 
the Pomeron. A feature of this model is the strong dependence on the wave functions 
of the involved hadrons, giving a possible source of uncertainty. The LN model is 
used to avoid the infrared divergences employing a modified gluon propagator, with 
the forms displayed above. 

The scattering amplitude can be written aJ^3 



A pp (s,t) = Sisa 2 s {T1-T2), 



(9) 
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where 



T x = J d 2 k V (| + k) V (| - k) G 2 p (q, 0), (10a) 

T 2 = J d 2 k V (| + k) V (| - k) G p («, ft - |) [2G p (g, 0) - G p (q, k - |)}ob) 

where s,t — —q 2 are the Mandelstam variables and G p (q, k) is a convolution of the 
proton wave functions, 



G p (q, k) = J d pda ip*[a,p) tp(a,p — k — aq), 

and can be related with the Dirac form factor as follows, 

G p (q,0) ^F^q 2 ) (12a) 

q 2 



(**-§)-« K 



(12b) 



where 

4m^ - 2.7!!/ 

\2 ' 



J~i(*)= ,2 . (!3) 



(4m2 - t) (1 - t/0.71) 

where m p is the proton mass. 

The terms in eqs. (|10all0b|) have a simple interpretation: the former comes from 
Feynman diagrams where the gluons are connected to the same quark in the proton. 
The last one com es from diagrams where the gluons are connected to different quarks 
in the protorP^. 

The total cross section is given by the optical theorem, of ot = *42 P ( S > 0)A S > as 
well as the differential cross section da jdt = \A2 P (s, t)\ 2 /16irs 2 . For low momentum 
transfer, the elastic differential cross section is fitted by an exponential expression, 
da/dt = Ae Bt = af ot e Bt /(l6ir), where B is the logarithmic slope, given by B = 
d/dt[\n(da/dt)}\ t= o- 

With the scattering amplitude above, a finite result for t he tot al cross section 
is obtained due to the cancellation of the infrared divergence eP^ I Otherwise, the 
differential cross section has an unphysical result (B(t — 0) — > oo) if the perturbative 
gluon propagator is employed. When the Landshoff-Nachtmann model is used, the 
contribution from the term 7~2 is negligible and the remaining infrared divergences 
are regularized by the modified gluon propagators. 

However, the two gluon description of the scattering amplitude is a crude approx- 
imation, which does not yield any energy dependence in the above cross sections. 
The reason is the absence of the gluon ladders in the s channel, as in the BFKL 
framework, which gives an energy dependence in the cross sections. To restore this 
dependence, we assume that the amplitude above only describes the energy inde- 
pendent part of the cross section and an extra term, introduced ad hoc gives the 
Regge energy behavior, namely, 

/ „ \ «p(*)-l 
A pp (s,t)^(-\ A^(s,t), 
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where ctp(t) — a(0) + c/(i) is the (soft) Pomeron trajectory (witfPl a(0) = 1.08 
and a'(t) ~ 0.25 GeV~ 2 ) and s is an energy scale. The extra Regge term gives the 
following cross sections, 







a P (0)-l 


Ctot = 




} °tot! 




V s 0, 


da 




2 ap (t)-2 d(j0 


~dt = 


V s o, 


i ~dT 



(15a) 



(15b) 

The above model was used by Refs. 171 1221 1241 where the non-perturbative propa- 
gator comes from different methods, discussed in the next s essi on. The experimental 
data are ISR results for elastic proton-proton scattering aJ251 = 53 GeV. 

In the previous model for the scattering amplitude, the running coupling con- 
stant was considered fixed. In Ref. [2EI the value of the coupling was determined 
thought the effective Pomeron coupling with the hadrons, which depends on the 
gluon propagator and has well defined experimental value. 

Motivated by the good results of our previous calculations with a gluon propaga- 
tor and a frozen running coupling^^, we modify the former model of pp scattering, 
considering now a frozen running coupling constant and a modified gluon propa- 
gator. We emphasize that this is a phenomenological modification, since a running 
coupling is a next leading order contribution and the model employed is a first or- 
der contribution to the Pomeron exchange. As usual in previous works, the scale 
of the running coupling constant is chosen as the incoming momentum into the 
quark-gluon vertex. 

After the above considerations, the scattering amplitude in the modified model 
then reads as 

AW(s, t) = S,s (£) QP(t) " / J* a s (§ + k) V (f + k) a s (§ - k) V (f - k 

G p (g,0)-G p (<z,fc-!)] 2 . (16) 

where the terms from eqs. ()10all and (|10b|l were combined as only one term and the 
running coupling terms are included in the integrand. 

The production of massive vector mesons is a process in which the Pomeron 
assumes an important role^- 1 - 1 , being a powerful laboratory for the study of the 
diffraction. The relevant process is 

7 W +p^V +p 

where V is the vector meson, p is the proton and 7 is the photon, being real (photo- 
production) or virtual (lepto-production) . The process can be described in the high 
energy regime in a factorized way in three steps: first, the photon splits into a pair 
of quark and anti-quark; second, the pair interacts with the proton through the 
exchange of a color singlet, the Pomeron (the meson has the same quantum numbers 
of the photon) and last, the anti-quark-quark pair recombines into a bound state, 
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Fig. 2. Contributing Feynman diagrams to the meson vector production by two gluons exchange. 



the vector meson. The scattering amplitude can be factorized in these three steps, 
as will be discussed below. Relevant kinematic variables are the photon virtuality, 
Q 2 , the meson mass, My, the square of the exchanged momentum, t, and the 
center-of-mass energy squared of 7p, W 2 . 

In the following, we employ the model of Cudell and Roye in the description 
of the production of p mesons in HI and ZEUS, with the same modifications already 
used in the pp elastic scattering, meaning the same set of modified gluon propagators 
and running coupling constant. 

A factorization scheme is employed, where three steps are considered, one of each 
with the following assumptions: (a) the interaction of the Pomeron with the proton, 
where the former interacts only with one valence quark of the proton, considering 
the Dirac elastic form factor for the proton; (b) the coupling of the Pomeron with 
the pair and the transition qq — > V, modeled assuming that the meson is in the 
lowest Fock state (two quarks) and without any Fermi momentum, allowing the use 
of a non-relativistic model for the qqV vertex. However, Cudell's models does not 
give any energy dependence for the cross sections and, for this reason, a Regge-like 
energy term is added ad hoc, due to the exchange of only two gluons. In this case, 
the (b) step above includes the possibility that the Pomeron (as two gluons) can 
be attached in different or in the same quark in the virtual pair, as shown in fig. El 
However, contributions from diagrams where the gluons of the Pomeron are attached 
with different quarks of the proton are suppressed in order to reproduce the quark- 
counting rule in the LN modeH An improvement considering that the interaction 
of the two gluons with the proton comes from two Feynman diagrams, one in which 
the gluons are attached with the same quark in the proton and another in which 
the gluons are attached to different quarks was given in Ref. 28] The contributions 
of these diagrams eliminate the infrared divergences as pointed out by Ref. 023 The 
last contribution, the only one from attachment to different quarks, is modeled by 
a modification of the Dirac elastic form factor, as discussed below. 

In this paper, we consider the simple model of Cudell and RoyerP^, with the 
same modifications used in the description of the pp scattering, supported by several 
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reasons. First, its simple form provides an easy and direct modification of the gluon 
propagator and running coupling constant. Thus, a modified gluon propagator and 
a frozen running coupling constant as well as a Regge term to res tore the e nergy 
dependence are employed, supported by a series of previous workJZHESEEl This 
formulation has the goal to test the modified gluon propagators in association with 
a frozen running coupling constant, not present in any previous model to analyze 
this process. 

A discussion about the terms that contribute to the unmodified scattering ampli- 
tude can be founded in the original worlPSI. The amplitude for the possible allowed 
helicities (transverse to transverse and longitudinal to longitudinal, the other cases 
being suppressed by w~ 2 ) is given by 



2 64 9clm™VVWjv 



A TT = iw 2 2 * ~ (18) 



1 



„2 M 



d 2 k [Fift)-g 2 (fc,fc-A)][k-(k-^ 
k 2 (k — A) 2 t - m 2 v - Q 2 + 4k • (k - A) 



-v 

All = -^——Att, (19) 
my 

where w 2 is the center-of-mass energy squared of the quark-photon system, g^ lra = 
£V47TO!cim is the electromagnetic coupling of different vector mesons: £ = l/v2 for 
the p meson; fy = 0.025my is the vector meson decay constant squared, with my 
the meson mass; Fi(t) is the elastic Dirac form factor, 

4m 2 - 2.79i 
Fl{t ' = (477i 2 -t)(l-</0.71) 2 ' 

with m p the proton mass and £ 2 (fc, k — A) = i*\(k 2 + (k — A) 2 - k • (k - A)). The 
last term comes from the diagrams that couples with different quarks inside the 
proton and, in the original work, avoi ds the infrared divergences. 
The differential cross section i J2SI 

da d<7T dcjL R r ■ ,9 , 9i 

where e is the polarization of the photon beam (e w 1 at HERA) and R is a factor, 
needed to restore the energy dependence, since two gluon exchange does not give 
such dependence in the cross section. 
The amplitude, eq. (JT5J, now reads 



A TT = iw —= „ — — — Id ka s (k)a s (k - A) 



,.2 64 g^myy/myfy ( s \ ap{t) 1 y ;2| 
V6 m v + Q 2 -t \s 0/ 



A' LL — A' TT . (23) 

my 
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Fig. 3. Results of the modified model, fitting to the pp elastic scattering data with y/s = 53 GeV 
(right) and yfs = 1.8 TeV (left) with different gluon propagators and frozen coupling constants. 
The parameters of the fits are given in the text. 



The expression <|22[1 can be rewritten in a more compact way, 



4' 



TT 



1 r( -a P (i)~\ 



f (A, Q ) , 



where ( = s/s Q with s = Q 2 + my-t, C = (64/V6)g^ m my V m vfv and / (A,Q 2 ) 
covers the remaining terms. The differential cross section comes from the optical 
theorem, 



da 



1 c3f 2(ap(t)-l) /2 ( A g 2) 

at 107T 



as well as, the total cross section, 

/ da 1 , 



Q 2 



J dt f{<*r{t)-l)f (A,Q 2 ). 



(25) 



(26) 



4. Results 

The differential cross section for pp scattering is calculated using the eq. (|16|l into 
eq. (|15b(l . varying the propagators as well as the running coupling constants, as 
pointed out above. The results of the calculation are displayed in the fig. EH for two 
set s of experimental data for distinct ce nter of mass energies: y/s — 53 GeV from 
ISEpS and y/s=1.8 TeV from TevatrorM 

Some comments over the results are addressed in the following. In the case of 
the pp scattering, we see the significant variation in the cross section when different 
propagators are employed. The slope of the curves are in accordance with data, due 
the ad hoc term in the amplitude. The energy independent part of the amplitude, 
modeled by the eq. © gives a contribution on the t behavior, on the intercept 
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Fig. 4. (a) Total cross section obtained with the eq. I2ril for p meson. (b)p photo-production 
differential cross section, compared with the HI and ZEUS data. 



(da I dt\t—o) of the curve and the normalization of the result. If we allow the vari- 
ation of the propagator parameters, is easy to fit the curves to the data, for both 
energies. The calculation is sensitive to the value of the parameters, giving larger 
modifications when the massive parameter is changed. In opposition, when the £ 
parameter is changed, the change in the result is significantly smaller than with the 
mass parameter. Thus, we will use the value of the £ fixed as 4. 

In the case of the Cornwall propagator, the propagator's parameter best value is 
nig = 0.53 GeV, which lies in the range of validity of the approximation employed 
(see above), and for the Habel/Gribov propagator, the value for the parameter is 
b = 0.3 GeV, close to the previous one, used in Ref. |S1 In all the other propaga- 
tors, the parameters are the original ones. In two cases, the Gorbar/Natale one and 
Aguilar/Natale one, we obtain that the original choice of the parameters underes- 
timate the results. When employing the Alkofer et al. propagator, eq. JSJl, and the 
corresponding running coupling constant, eq. JJJ), we obtain a result that overesti- 
mates the experimental data for both energies. Similar result was found by Refs.EI 
0for another observables, supporting our results. 

The results of the modified model for the total cross section, eq. 1)26(1 . are dis- 
played in figure 01 for the gluon propagators and coupling constants discu ssed in 
the previous section and in comparison with the HERA (HI and ZEUS) dat tP^ffi 

The center-of-mass energy is w ~ 100 GeV, the result being less sensitive to this 
parameter. Due to the reasons exposed above, the energy dependent term uses the 
soft Pomeron trajectory, ap(t) = a(0) + a'(t), with a(0) = 1.08 and a'(t) ~ 0.25 
GeV -2 . As in the pp elastic scattering case, the £ parameter of the running coupling 
constant is fixed in £ = 4, its change does not modify significantly the results. 

Therefore, the gluon parameters are, for the cases with some freedom of choice: 
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Cornwall, m g — 0.4 GeV and Habel/Gribov, b = 0.2 GeV. The final result is very 
sensitive to this parameter, giving very different results if altered by a small amount, 
allowing the fit of the data, if the parameters are free. 

The results lie along in the same general features of the pp scattering: if the 
original values of the propagator parameters are employed, the data are not fitted, 
with the exception of two propagators, which have some choice of the parameter. 
The overall behavior describes the data, although the normalization which depends 
on the choice of the propagator is large. However, if the propagator parameter is 
allowed to vary, we can improve the agreement of the fit, as in the pp scattering, but 
in the case of the total cross section, in the high-Q 2 region, the result stays above 
the data, showing that the model can describe better the kinematic region where 
the non-perturbative effects are relevant. The differential cross section of photo- 
production (Q 2 — 0) of the p meson is displayed in the fig. 01 where the differences 
between the results from different propagators are evident, but we have a reasonable 
description of the data. As in the proton-proton case, some of the propagators lie 
below the data and the Alkofer's choice lies above the cross section data. 

5. Conclusions 

We present a comparison of different results for the gluon propagator and the related 
frozen running coupling constant, which functional forms are obtained from different 
non-perturbative methods. These results are tested in an appropriated framework, 
the high energy diffractive scattering, where the low momentum transfer effects are 
relevant, and therefore the application of the modified gluon propagator is relevant. 
Among the several processes, our choices are the elastic pp scattering and light me- 
son production. The same ideas were employed by the present authors, in different 
processes also with a sma ll momentum transfer as, for example, the description 
of meson form factorJ^E] anc j the process 77 — > J/t[)J/t[) also was described in a 
similar framework in Ref. |H| 

The original idea is not to perform a fit to the data, but confront the original 
solutions for gluon propagators through specific observables. The different values for 
the massive parameter obtained for the pp elastic scattering process indicates some 
weakness in this approach. The same problems appear in the meson production, 
where the massive parameter is different of the pp case. If the functional forms for 
the gluon propagator is correct, what is the origin of these discrepancies? An answer 
could be the model employed in the calculation, which have the same common origin, 
but then the results would be similar, which is not the case. 

The results for the description of the observables are not quite satisfactory, 
showing that the approximation employed here, although still crude, signs that 
the use of the modified gluon propagator in diffractive events can be a reasonable 
description of the data, although not complete. The crudeness of the fit in some cases 
can indicate an incomplete picture for diffractive scattering, and also our choice of 
not performing a fit to the experimental data, adjusting the massive parameter, may 
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be a very restrictive one. Moreover, the present results show that the assumptions 
of the modified model should be used with caution. The gluon propagators are 
obtained in a self consistent way, although the original choices for the parameters do 
not give good results. Otherwise, the model employed can be incomplete, omitting 
contributions from other Green functions, as quark propagators, for example. Thus, 
a possible natural extension is to employ a complete set of Green functions: the 
quark and gluon propagators and the quark-gluon vertex. 

In the present processes, in which the choices have the feature of a dynamical 
mass (or a massive parameter), we obtain a quite good agreement with experimen- 
tal data for distinct energies, for two choices for propagators, the Cornwall's and 
Habel/Gribov ones, which have some freedom to vary the parameters. The results 
obtained are along the same lines of the previous attempts to describe these pro- 
cesses, however considering distinct propagators. The results show that the employ 
of modified gluon propagators is quite valuable in phenomenology, although requires 
some caution, since the model used to construct the scattering amplitude is a first 
approximation. The results point out to a gluon propagator with dynamical mass 
and frozen running coupling constant, although the present results in addition with 
the previous one J^I'^ m do not definitely determine the most suitable expression 
for the propagator, mainly due the approximations employed and the still poor 
knowledge on the interplay of the whole aspects involved in these calculations. 
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